arXiv:1503.03656vl [hep-lat] 12 Mar 2015 
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We present the results of a lattice study of the second moment of the light-cone pion distribution amplitude 
using two flavors of dynamical (clover) fermions on lattices of different volumes and pion masses down to 
m,r ~ 150 MeV. At lattice spacings between 0.06 fm and 0.08 fm we find for the second Gegenbauer moment 
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I. INTRODUCTION 

Hard exclusive processes involving energetic pions in the 
final state are sensitive to the momentum fraction distribution 
of the valence quarks at small transverse separations, usually 
called the pion distribution amplitude (DA). Classical appli¬ 
cations have been to exclusive two-photon processes, 

e.g., the pion electromagnetic form factor at large momentum 
transfer and the transition form factor 7 * ttj. The latter 
process plays a very special role as the simplest hard exclusive 
reaction where QCD factorization can be tested at a quanti¬ 
tative level. It received a lot of interest recently, triggered 
by the partially conflicting measurements by BaBar ID and 
BELLE 13 Up to photon virtualities of the order of 40 GeV^, 
see, e.g., isHini. Arguably, the most important application of 
the pion DA is currently the study of semileptonic weak de¬ 
cays B —r irfDe at large recoil lfTTI4T3]l using light-cone sum 
rules (LCSR) HI US) and weak hadronic decays B -4- tttt 
etc. in the framework of QCD factorization ifT^ fTTll . Both 
reactions contribute prominently to the determination of pa¬ 
rameters of the quark mixing matrix in the Standard Model. 

The precise definition of the pion DA is based 

on the representation GHS as the matrix element of a non¬ 
local light-ray quark-antiquark operator. Lor example, for a 
positively charged pion 

{0\d{z2n)i^j5[z2n, zin]u(zin)|7r(p)) 

= tfAp-n) (1) 

Jo 

where is the pion momentum, is a light-like vec¬ 
tor, v? = 0, Zi ^2 are real numbers, [z 2 n, Zin] is the Wil¬ 
son line connecting the quark and the antiquark fields and 
= 132 MeV is the usual pion decay constant. The DA 
4>Tr{x, is scale-dependent, which is indicated by the argu¬ 
ment p?. 

The physical interpretation of the variable x is that the u- 
quark carries the fraction x of the pion momentum, so that 
1 — X is the momentum fraction carried by the d-antiquark. 


* Paula.Perez-Rubio@ur.de 


Neglecting isospin breaking effects and electromagnetic cor¬ 
rections the pion DA is symmetric under the interchange 

X 44 1 — x: 


4>^{x,p?) = x,p'^). ( 2 ) 

Due to this symmetry, only the even moments involving the 
momentum fraction difference 

^ = X — (1 — x) = 2x — 1 (3) 


carry nontrivial physical information: 


iO = [ dx{2x-lff(j)^{x,p^), n = 0,2,.... (4) 

Jo 

The definition in Q implies the normalization condition 



( 5 ) 


A convenient parameterization of DAs is provided by the 
conformal expansion llT8l420l . The underlying idea is to use 
the conformal symmetry of the QCD Lagrangian to separate 
transverse and longitudinal variables in the light-front pion 
wave function, similar in spirit to the partial-wave decom¬ 
position in quantum mechanics. The dependence on trans¬ 
verse coordinates is formulated as a scale dependence of the 
relevant operators and is governed by renormalization-group 
equations. The dependence on the longitudinal momentum 
fractions is described in terms of Gegenbauer polynomials 
C'^^^( 2 x — 1 ) which correspond to irreducible representations 
of the collinear conformal group SL(2,IR). In this way one 
obtains 


p^) = 6 x(l — x) 


1 + ^ affp^)Cl/\2x-l) 


2 = 2 , 4 ., 


( 6 ) 

where all nonperturbative information is contained in the set 
of coefficients (Gegenbauer moments) a„(/j,g) at a certain 
reference scale pQ. To leading-logarithmic accuracy (LO), 
the Gegenbauer moments renormalize multiplicatively with 
the anomalous dimensions rising slowly with n. Thus the 
higher-order contributions in the Gegenbauer expansion are 
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suppressed at large scales and asymptotically only the leading 
term survives. 


= 6 x(l - x), (7) 

which is usually referred to as the asymptotic pion DA. It is 
widely accepted, however, that the pion DA deviates signifi¬ 
cantly from its asymptotic form at scales that can be achieved 
in experiments. 

A particular model of the pion DA proposed by Chernyak 
and Zhitnitsky in 1982 121] has played an important role in 
historic perspective. It was based on a calculation of using 
QCD sum rules fT2\ . which resulted in a large value 02 
0.5 — 0.6 (at the scale 1 GeV), and the assumption that all 
higher-order coefficients can be neglected. 

Since then, different approaches have been used; QCD sum 
rales with various improvements (e.g. ||23H23), LCSR-based 
analysis of experimental data on the pion electromagnetic and 
transition form factors (e.g. C®) and weak S-meson de¬ 
cay form factors (e.g. IH), lattice calculations EH Ell and 
recently also in the framework of Dyson-Schwinger equa¬ 
tions EH- A recent compilation of the existing results for 
02 can be found in Table 1 of Ref. Q. 

Estimates of yet higher-order Gegenbauer coefficients are 
rather uncertain. A direct calculation of 04 proves to be dif¬ 
ficult and its extraction from the experimental data on, e.g., 
the pion transition form factor is complicated by the fact the 
LO contribution is proportional to the sum of Gegenbauer mo¬ 
ments 

/ —= 3[1 + + a4(p^) + ■ ■ ■] ■ (8) 

Jo X 

Thus, the values of a 2 {n'^) and a 4 (/i^) obtained in these ex¬ 
tractions appear to be strongly correlated. The strong scal¬ 
ing violation in the pion transition form factor observed by 
BaBar III (but not confirmed by BELLE Q) would imply a 
considerable enhancement of the pion DA close to the end¬ 
points, meaning that the expansion in Gegenbauer polynomi¬ 
als is converging very slowly if at all, see the detailed discus¬ 
sion in Eiiiiia. The forthcoming upgrade of the Belle exper¬ 
iment and the KEKB accelerator E9l , which aims to increase 
the experimental data set by a factor of 50, will allow one to 
measure transition form factors and related observables with 
unprecedented precision and resolve this issue. The question 
at stake is whether hard exclusive hadronic reactions are under 
theoretical control, which is highly relevant for all future high- 
intensity, medium energy experiments like, e.g., PANDA. On 
the theory side, several proposals exist how it might be pos¬ 
sible to access DA moments beyond the second one (or the 
DA pointwise in x) on the lattice, e.g., MM, but the corre¬ 
sponding techniques are only in the exploratory stage. 

In this work we extend the lattice study lf26l of the sec¬ 
ond moment of the pion DA by making use of a larger set 
of lattices with different volumes, lattice spacings and pion 
masses down to ^ 150 MeV and implementing sev¬ 
eral technical improvements. We employ the variational ap¬ 
proach with two and three interpolators to improve the signal 
from the pion state. The renormalization of the lattice data 


is performed nonperturbatively utilizing a version of the RT- 
SMOM scheme. Eor the first time we include a nonpertur- 
bative calculation of the renormalization factor corresponding 
to the mixing with total derivatives, which proves to have a 
significant effect. Our main result is 

02 = 0.1364(154)(145)(?) (9) 


for the second Gegenbauer moment of the pion DA, and 


= 0.2361(41)(39)(?). 


( 10 ) 


Both numbers refer to the scale p, = 2 GeV in the MS 
scheme. The first error combines the statistical uncertainty 
and the uncertainty of the chiral extrapolation. The second er¬ 
ror is the estimated uncertainty contributed by the nonpertur- 
bative determination of the renormalization and mixing fac¬ 
tors. Our lattice data are collected for the lattice spacing 
a = 0.06 — 0.08 fm, and this range is not large enough to 
ensure a reliable continuum extrapolation. The correspond¬ 
ing remaining uncertainty is indicated as (?). It has to be ad¬ 
dressed in a future study. 

The paper is organized as follows. In the next section we 
discuss the aspects of the continuum description of the pion 
DA that are relevant for our work. The basics of the lattice for¬ 


mulation are given in Sec. Ill An important ingredient in our 
calculation is the nonperturbative evaluation of the renormal¬ 
ization and mixing coefficients, which is described in Sec.|IV| 
The methods applied in the analysis of the bare data are de¬ 
tailed in Sec.|V] Our results are presented in Sec.|VI[ followed 
by our conclusions and an outlook. In an Appendix we collect 
Tables of intermediate results for each gauge field ensemble 
used in our work. 


II. MOMENTS OF THE PION DISTRIBUTION 
AMPLITUDE 


The nonlocal operator in the expression for the pion DA ([T]) 
is defined as a generating function for renormalized leading- 
twist (i.e., twist two) local operators. 


d{z2n)i^"f5[z2n, zin]u{zin) = 
/c I 

^ kill 

k,l—0 


( 11 ) 


where 


= ... 79^;,£>^,+1 ... 79 ^,+, 7 p ) 75 m ( 0 ) . ( 12 ) 

Here 79^ is the covariant derivative and (...) denotes the sym- 
metrization of all enclosed Lorentz indices and the subtraction 

(k 1) 

of traces. The local operators M.pdi...fi.k+i ^re assumed to be 
renormalized, e.g., in the MS scheme. 

As a consequence, moments of the pion DA are given by 
matrix elements of local operators: 


= */7rP(pPMi • ..Pp,^,){x\l-X)'^). 

( 13 ) 
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Neglecting isospin breaking effects and electromagnetic cor¬ 
rections one obtains the symmetry relation 

( 14 ) 

and thus 

{x\l-x)’^) = {x'^{l-xy). (15) 


In addition, the product (Leibniz) rule for derivatives 

■I + I + i ^-I ‘ppi...pfc + i + i + ‘ppi.,,ps,_,.,) 

gives rise to the momentum-conservation constraint 

(a;'+i(l-tc)'=) + (a;'(l-x)'=+i) = (x'(l-a;)'=) . (17) 


(16) 


Specializing to the second moment, Z -|- fc = 2, it is easy to see 
that only one independent matrix element remains, e.g., 

(^ 2 ) = 1_4(x(1-x)) (18) 


or 

a 2 = ^(C|/"( 2 x-l)) = ^[5(a-l] 

= ^[l-5(x(l-x))] , (19) 

so that any moment (^^), 02 , (x(l — x)) etc. can be used as 
a nonperturbative parameter to characterize the shape of the 
pion DA. Lacking any a priori information on the relative size 
of the different contributions, all such choices are equivalent. 
It is widely expected, however, that the numerical value of 
is not far from 1/5 corresponding to the asymptotic pion 
DA 0. Hence, if 

( 20 ) 

is determined with a given accuracy at some reference scale 
/io by a certain nonperturbative method, and 02 is then ob¬ 
tained from the relation ( | 20 ] l, the error on is strongly am¬ 
plified by the subtraction of the asymptotic contribution. This 
effect is well known and has been observed in all calculations 
up to date. The error on 02 is relevant as it propagates through 
the renormalization group equations. In other words, although 
using 02 as a nonperturbative parameter instead of for the 
pion DA at a low reference scale (pT^ (x, /ig) is just a rewriting, 
this choice is much more adequate in order to describe the 
pion DA at high scales, (pT^{x^Q‘^), Q ^ p-o, which enters 
QCD factorization theorems. Another issue to consider is that 
the relation in Eq. and therefore ( [T 8 | , ( [T^ , ( [ 20 | l can be 
broken by lattice artifacts. Thus the choice of suitable opera¬ 
tors requires some care. We will discuss our choice in more 
detail in the next section. 


III. LATTICE FORMULATION 


define the bare operators 


= d{x) 

= d{x) 


D{pDu — 2D(^pD^ + 

D{p,Dy + 2D(^pD^ + D(^pDy 


7p)75 u{x ), 

7p)75'u(x) 

( 21 ) 


as our operator basis. On the lattice the covariant derivatives 
will be replaced by their discretized versions. 

The operator 0~^^, can be written in a conventional short¬ 
hand notation as 


= d{x)D(^pD^-ip)-f^ u{x) (22) 

and its matrix element between the vacuum and the pion state 
is proportional to the bare lattice value of ((x — (1 — x))"^) = 
(^^). In the continuum, the operator is the second 

derivative of the axial-vector current; 

Ojp^ix) = di^pd^Op){x) with Op{x) = d{x)jp-f 5 u{x). 

(23) 

However, this relation is violated on the lattice because of 
discretization errors in the derivatives. The distinction be¬ 
tween and di^pd^Op) for finite lattice spacing appears 

to be numerically important and will be discussed in detail in 
what follows. Note that C9+ is the Euclidean analogue of the 
Minkowski-space operator Ad^uch 
that its matrix element between the vacuum and the pion state 
corresponds to the bare value of ((x -f 1 — x)^) = ( 1 ^). 

The corresponding renormalized (e.g., in the MS scheme) 
axial-vector current is then given by 

0f^{x) = ZaOp(x) (24) 

with Za ^ 1 on the lattice. 

In order to express its matrix elements in terms of the phys¬ 
ical quantities introduced in Minkowski space we apply the 
rules 


7m = 74 , 7m = *7i (25) 

for j = 1,2,3, where the subscript M distinguishes the 
Minkowski objects. Consequently, 

7“ = *7m7m7m7m = -71727374 = -75 ■ (26) 

The components of the three-vector p of the spatial momen¬ 
tum of the pion will be denoted by pj , although they are equal 
to the contravariant space components of the Minkowski mo¬ 
mentum p. The time component of the Minkowski momentum 
is identified with the corresponding energy: po = i7,r(p)- In 
this way one gets in Euclidean notation 

(0|Or(0)k(p)) = -iE^ip)U , (27) 

(0|OfS(0)|7r(p)) = -p,^. (28) 


While the above relations refer to renormalized operators 
in Minkowski space, we now move to Euclidean space and 


Similarly, the Euclidean space components of the coordi¬ 
nate vector X are identified with the contravariant components 
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of the Minkowski space-time four-vector, while for the time 
components we have xq = —ix^. This entails the following 
rule for the covariant derivatives: 

- iDf = Di , Df = Dj . (29) 

Therefore we find, e.g., for j ^ k 

(0|C>“®"(0)|7r(p)) = iU{^'^)E^{p)pjPk . (30) 

The operators and mix under renormalization 

even in the continuum. On the lattice the continuous rota¬ 
tional 0(4) symmetry of Euclidean space is broken and re¬ 
duced to the discrete H (4) symmetry of the hypercubic lattice. 
This symmetry breaking can introduce additional mixing op¬ 
erators. It can even lead to mixing of the operators of interest 
with operators of lower dimension such that the mixing coeffi¬ 
cients are proportional to powers of 1 /a. This complicates the 
renormalization procedure significantly. However, it may be 
possible to choose the lattice operators such that they belong 
to an irreducible representation of H (4) which forbids mixing 
with further operators, in particular with lower-dimensional 
operators. In the present case there is one such choice, given 
by the operators with all three indices different. For the 
computation of the required matrix elements we can restrict 
ourselves to the operators (see, e.g., mwi) 

^^4 %. jV fee {1,2,3}. (31) 

The renormalized operators are then given by 

(^) + Z, 20 t,M , 

of+{x) = Z^20t,,{x). (32) 


Note that due to the discretization artifacts in the derivatives 
one cannot expect Z 22 to be equal to Za- 

For the calculation of q,ms .^g ^j.g jg^j 

two tasks: computation of the bare matrix elements and eval¬ 
uation of the renormalization factors. We extract the bare ma¬ 
trix elements from two-point correlation functions of the oper¬ 
ators and Op with suitable interpolating fields J{x) for 
the TT-mesons. For the latter we consider the two possibilities 

Jz{x) = u{x)"f^d{x ), 

Ji^{x) =u{x)"iA'y^d{x) (33) 

with smeared quark fields. The details of our smearing algo¬ 
rithm will be given below. Fet 


p) = ^ e-*P "(0,(x, f) J^(0)), 

X 

,(x,f)J^(0)), (34) 


where A = 5 01 A = 45, p is the three-vector of the spatial 
momentum, and the summation goes over the set of spatial 
lattice points x for a given Euclidean time t. 

For times t, where the correlation functions are saturated by 
the contribution of the lowest-mass pion state, we expect that, 

e.g.. 


= (0|C)± (0)|7r(p))(7r(p)|J^(0)|0) 






2E I 




(35) 


Here E = F^,n.(p), T is the temporal extent of our lattice, 
and the r-factors take into account transformation properties 
of the correlation functions under time reversal. One finds 
Tjs = -1, = 1 , to = 1 for the operators O 4 

and To = — 1 for Oj, where j,k = 1,2,3. We utilize these 
symmetries in order to reduce the statistical fluctuations of our 
raw data, i.e., we average over the two corresponding times t 
and T — t with the appropriate sign factors. 

From the ratios 


±.A _ G^i’^{t,p) 


(36) 


we can extract the required bare matrix elements 
(0|(!lp^^(0)|7r(p)), which carry the information on the 
second moment of the pion DA. 

Equation ( (30l l shows that a calculation of matrix elements 
of O^jf, requires two nonvanishing spatial components of the 


momentum. We choose them as small as possible, p = 27r/L, 
where L is the spatial extent of our lattice. To suppress sta¬ 
tistical fluctuations we average over the possible directions, 
e.g., p = (p,p,0), p = (p,-p,0), p = (-p,p, 0 ), p = 
(—p, —p, 0) for j = 1, fc = 2. If the correlation functions are 
dominated by the single-pion states, the time-dependent fac¬ 
tors in the ratios of correlation functions cancel and we obtain, 
e.g., for the operator 0 ^i 2 the momentum p = (p,p, 0) 

^«214 = -(y) 

where the constants are related to the bare lattice values 
of the second moment of the pion DA through 

(^2)b“e = flbare ^ J_ ^5^- _ j^+'^ (33^ 

They should not depend on the choice of the interpolating held 
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Ja- Note that and therefore for bare quantities 


^bare 


^_(5(^2)bare-i) ^ 


(39) 


For the renormalized moments in the MS scheme we obtain 




^MS 


12 


5Cll^ + (5 Ci2 — C22)^~'' 


(40) 


where 


> _ -^11 /■ _ ^12 /- _ Z22 

bit — J (,12 — 5 (,22 — 

ZjA 2 :^ 


(41) 


are ratios of renormalization constants defined in the next sec¬ 
tion. 

In the continuum limit we expect that 

2'22(0|C9+fc(0)|7r(p)) = -ZAPjPk{0\Oi{{))\'n{v)) 

= ipjPkE^{p)f7v ■ (42) 


Hence the quantity 




1/22 


(0|O+,(0)|7r(p)) 


(-PjFfe)(0|(94(0)|7r(p)) 


= C22i?+ (43) 


should approach unity as the lattice spacing tends to zero. In 
this case the relation 


a. 


™ = 2(5(a“"-i) 


(44) 


is recovered (cf. Eq. ®), whereas for finite lattice spacing it 
follows from ( |40| l 


ar = ■ 


(45) 


We emphasize that Eq. (44 1 is only recovered in the continuum 
limit, which is always delicate. There are two possibilities; 
Either is measured on the lattice, the result extrapolated 
to zero lattice spacing, and at the final step 02 is obtained using 
the relation ( |44| i, or 02 is calculated directly on the lattice and 
then extrapolated to the continuum limit. The first approach 
was used in Refs. Il26l|22l whereas in this work we use the 
second method. 


IV. RENORMALIZATION CONSTANTS 


Erom our bare matrix elements we have to compute the cor¬ 
responding renormalized matrix elements in the MS scheme, 
which is used in the perturbative calculations. In the contin¬ 
uum we therefore have to deal with the renormalization of the 
two mixing operator multiplets given in Eq. 0. Note that 
(^p)ii/, being the second derivative of the axial-vector current, 
has vanishing forward matrix elements, at least in the contin¬ 
uum. 


TABLE I. Ensembles used for nonperturbative renormalization. 


p 

K 

Size 

5.20 

0.13550 

32® X 64 

5.20 

0.13584 

32® X 64 

5.20 

0.13596 

32® X 64 

5.29 

0.13620 

32® X 64 

5.29 

0.13632 

32® X 64 

5.29 

0.13640 

64® X 64 

5.40 

0.13640 

32® X 64 

5.40 

0.13647 

32® X 64 

5.40 

0.13660 

48® X 64 


On the lattice we work with the operator multiplets 

^^23 ’ ^413 I ^412 J ^123 (46) 

and 

6^423 ’ ^il3 > ^412 > 69 i 23 • ( 47 ) 

Under the hypercubic group iT(4), both multiplets transform 
identically according to a four-dimensional irreducible repre¬ 
sentation ll3^ . The symmetry properties of these multiplets 
ensure that they do not mix with any other operators. Because 
of the well-known shortcomings of lattice perturbation theory 
we want to determine the renormalization and mixing factors 
nonperturbatively on the lattice, utilizing a variant of the RT- 
MOM scheme. However, since forward matrix elements of 
eventually vanish, we cannot use the momentum geom¬ 
etry of the original RI’-MOM scheme but have to work with a 
kind of RI’-SMOM scheme 0^ . 

In order to describe our renormalization procedure we con¬ 
sider a somewhat more general situation than what is needed 
in this paper. Let (x) (i = 1,2,... ,d,m = 1,2,..., M) 
denote M multiplets of local quark-antiquark operators which 
transform identically according to some irreducible, unitary, 
d-dimensional representation of id(4). Call the unrenormal¬ 
ized, but (lattice-)regularized vertex functions (in the Landau 
gauge) V^'^\p, q), where p and q are the external quark mo¬ 
menta. The corresponding renormalized (in the MS scheme) 
vertex functions are denoted by (7^™^ g). The dependence 

of on the renormalization scale p is suppressed for 

brevity. Note that 1^^™^ carries Dirac indices and is therefore 
to be considered as a 4 x 4-matrix. (The color indices have 
been averaged over.) 

We choose 


p= ^(1,1, 0,0) , 9= ^(0,1,1,0) (48) 

such that p^ = q^ = {p — qY = As our renormalization 













6 


TABLE II. Choices for the fits. 


Fit 

number 

Fit interval 
(in GeV^) 

^loops 

Lattice 

artifacts 

ro 

(in fm) 

roA^jg 

1 

4 < < 100 

2 

As 7^ 0 

0.50 

0.789 

2 

2 < < 30 

2 

As ^ 0 

0.50 

0.789 

3 

4 < < 100 

1 

As ^ 0 

0.50 

0.789 

4 

4 < < 100 

2 

As=0 

0.50 

0.789 

5 

4 < < 100 

2 

As 0 

0.49 

0.789 

6 

4 < < 100 

2 

As ^0 

0.50 

0.737 



TABLE III. Fit results at /3 = 

5.40 for /ig 

= 4GeVA 


Fit 1 

Fit 2 

Fit 3 

Fit 4 

Fit 5 Fit 6 

Cii 

2.026 

2.031 

2.123 

2.001 

2.040 2.041 

Cl 2 

-0.199 

-0.205 

-0.233 

-0.188 

-0.202 -0.203 

C22 

1.474 

1.476 

1.479 

1.467 

1.474 1.474 


condition we take (in the chiral limit) 

M d 

= Z-^ ^ Zw' E , (49) 

m"—l i—1 

where is the lattice Born term corresponding to . 
The wave function renormalization constant of the quark 
fields Zq is determined from the quark propagator, as 
usual ll?4ll . and subsequently converted to the MS scheme. Us¬ 
ing the lattice Born term instead of the continuum Born term 
and proceeding analogously in the calculation of Zq ensures 
that Z is the unit matrix in the free case. 

The renormalization matrix Z leads from the bare operators 
on the lattice to renormalized operators in our SMOM scheme. 
The matrix Z transforming the bare operators into renormal¬ 
ized operators in the MS scheme is then given hy Z = CZ, 
where the matrix C is defined as 

M d d 

m" — l i—1 i—1 

Here is the renormalized vertex function in the MS 

scheme and is the continuum Born term such that the 


conversion matrix C is completely determined from a contin¬ 
uum calculation. 

Here we have to consider the cases M = 2, d = 4 for the 
multiplets ( [46| ), ( [47| and M — 1, d = 4 for the axial-vector 
current. The required MS vertex functions in the chiral limit 
for up to two loops can be extracted from Refs. Il35l[3^ . As 
we are only interested in ratios of renormalization factors, Zq 
drops out and is not needed. In the following we describe our 
method for the determination of the renormalization matrix of 
the multiplets ( |46] ), ( [47] ). The procedure for the ratios with Za 
is completely analogous, because the anomalous dimension of 
the nonsinglet axial-vector current vanishes. 

The calculation of the vertex functions with the help of mo¬ 
mentum sources is straightforward. Partially twisted bound¬ 
ary conditions applied to the quark propagators allow us to 
vary the renormalization scale p independently of the lattice 
size. The ensembles used for the evaluation of the Z matrices 
according to the above formulas are listed in Table |I] Due to 
the rather small quark masses the subsequent chiral extrapo¬ 
lation appears to be quite safe. 

Ideally, the renormalization scale p should satisfy the con¬ 
ditions 


1/L^ < A|cd < < 1/a^ (51) 

for a lattice with lattice spacing a and extent L. Then lat¬ 
tice artifacts would be negligible and the scale dependence 
could be described by low-order continuum perturbation the¬ 
ory. However, the above conditions are hard to realize in prac¬ 
tice and the Z-values at any given scale suffer from discretiza¬ 
tion artifacts as well as from truncation errors of the perturba¬ 
tive expansions. Therefore we try to exploit as much of the 
available nonperturbative information as possible by perform¬ 
ing a joint fit of the /i-dependence of the chirally extrapolated 
renormalization matrices Z{a,fj,)MC for our three ^-values 
13 = 5.20, 5.29 and 5.40. 

The choice of the fitting procedure is motivated by the fol¬ 
lowing considerations. The (perturbative) running of the Z- 
matrices is governed by the anomalous dimension matrix 

^ = - ("f) ^ 

Introducing the running renormalized coupling p(/r) with 
^J.dg/d^x = P{g) we get 


^ = T(g) 7 

dg 13 (g) 


(53) 


This system of differential equations can formally be solved 
in the form 


z{fd)z-\^^o) = Y.i-^y 


raM 


n—0 


'g{ao) 


dgn 


's(mo) 


dgn- 


1 • 



7(g2) 7(gi) 

/3(g2) /3(gi) ’ 


(54) 


From the three-loop anomalous dimension matrix one can calculate a corresponding approximation of 
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FIG. 1. Renormalization and mixing factors (^ij in the chiral limit, perturbatively scaled to /tq = 2 GeV (cf. Eq. l |56[ l) together with curves 
representing Fit 1. The error bars show the statistical errors. The horizontal lines indicate the fitted values C,ij (a, /tq ). Note that the fit is aimed 
at describing the data for large values of the scale /r, the fit interval being 4 GeV^ < < 100 GeV^ 


which should describe the /i-dependence for 
sufficiently large scales /r, if there were no discretization ef¬ 
fects. Adding a plausible ansatz for an effective description of 
these lattice artifacts we arrive at the following fit function for 
the matrices Z[a, 

Z{a,n)Mc = VF(Af,Ato)^(a,/ro) -f 

+ A3(aV"f- (55) 

The fit parameters are the entries of the three renormalization 
matrices Z{a, fio) at the reference scale /xq and the entries of 
the three matrices Ai parameterizing the lattice artifacts. Note 
that we allow for a nonvanishing value of Z 21 although Z 21 
vanishes in the continuum. 

The statistical errors of the data are quite small, in particu¬ 
lar for larger scales, and the resulting statistical errors of the 
fit parameters turn out to be unrealistically tiny. Therefore 
the statistical errors will be ignored in the following. The 
systematic uncertainties, on the other hand, are much more 
important. In order to estimate them we perform a number 
of fits varying exactly one element of the analysis at a time. 
More precisely, we choose as representative examples for fit 


intervals 4 GeV^ < p? < 100 GeV^ and 2 GeV^ < fj? < 
30 GeV^, and we use the expressions for the MS vertex func¬ 
tions with nioops = lj2. For the parameterization of 

the lattice artifacts we either take the complete expression in 
Eq. ( [55] l or we set A 3 = 0. Finally, we consider values for 
tq and roAj^ corresponding to the results given in Ref. llJTl . 
The various possibilities are compiled in Table [n| 


As an example we show the fit results for /? = 5.40 in Ta- 
choosing /Xg = 4 GeV^. The numbers for the other 


III 


ble 


/3-values are similar. 


The largest effect comes from the variation of nioops: Work¬ 
ing with the 1 -loop vertex functions increases the result for 
(^11 by about 5%, and the modulus of the mixing coefficient 
(^12 increases even by about 17%. In order to obtain our final 
numbers for and extract them from 

the raw data for using each of these sets of values for Cn, 
C 12 and C. 22 - So we get six results for each of our gauge field 
ensembles. As our central values we take the results from Fit 
1. Defining Si as the difference between the result obtained 
with the <^s from Fit i and the result determined with the (s 
from Fit 1, we estimate the systematic uncertainties due to the 

















































8 


renormalization factors as •\/S 2 + (0.5 ■ + 6l + + 6 q. 

Here we have multiplied ^3 by 1 /2, because going from two 
loops to three or more loops in the perturbative vertex func¬ 
tions is expected to lead to a smaller change than going from 
one loop to two loops. This should amount to a rather conser¬ 
vative error estimate. 

In Fig. [^we show the entries of the matrix 


X 


AioV + 


(56) 


for Hq = 4GeV^ at our three /3-values along with the fit 
curves resulting from Fit 1 in Table [I^ The horizontal lines 
represent the fitted values Cii(a, Mo) etc. 

In the previous paper ll26l the renormalization and 
mixing factors were evaluated in a mixed perturbative- 
nonperturbative approach, based on the representation of 
0 + as the second derivative of the axial-vector current (see 
Eq. (p^). Repeating this calculation in a completely nonper- 
turbative setting we find that the overall renormalization fac¬ 
tor corresponding to Cii agrees within a few percent. The 
nonperturbative mixing coefficient, on the other hand, has the 
same (negative) sign as its perturbatively computed counter¬ 
part, but its modulus is up to one order of magnitude larger. 
This observation underlines the necessity of nonperturbative 
renormalization, at least for the presently reachable /3-values. 


V. ANALYSIS OF THE BARE DATA 


As was already mentioned in Sec. Ill the bare matrix el¬ 
ements related to the pion DA’s second moments can be ex¬ 
tracted from ratios of lattice correlation functions given by 
Eq. p 6 | ). We briefly describe our procedure. 

The gauge field configurations used in this work have been 
generated with the Wilson gauge action and Uf — 2 flavors 
of nonperturbatively improved Wilson fermions. We have 
analyzed 0(1000 — 2000 ) configurations for three different 
values of the gauge coupling, /3 = 5.20, 5.29, 5.40, and pion 
masses in the range m,r 500 — 150 MeV. The lattice spac- 
ings and spatial volumes vary between 0.06 — 0.081 fm and 
(1.71 — 4.57 fm)^, respectively. A list of our ensembles can 
be found in Table Eor scale setting we used the Sommer 
parameter with the value tq = 0.5 fm lEillMl- 

The correlation functions ( [34] l have been computed for 
the operators O 4 , O^jk (see Eq. ( [3T] i) leading to the ratios 
where Jopt is discussed below. On most of the en¬ 
sembles, we performed more than one measurement per con¬ 
figuration to increase the statistics. The source positions for 
the correlation functions were selected randomly to reduce the 
autocorrelations among configurations lying close to one an¬ 
other in the Monte Carlo history. We want the interpolating 
operators to have a good overlap with the ground state of the 
pion. To this end, Wuppertal smearing Il39ll40ll was applied to 
the sources, with APE smeared ED gauge fields. 

In order to reduce the overlap with excited states even fur¬ 
ther we have used the variational method ||42H45]| with the 


TABLE IV. Ensembles used for this work. 


n 

m-^l MeV 

Size 

m-n-L 

Number of 
configspi 


J = 5.20, a = 

0.081 fm, a~^ 

= 2400 MeV 


o.issge'^ 

280 

32® X 64 

3.7 

1999(x4) 


J = 5.29, a = 

0.071 fm,a“® 

= 2800 MeV 


0.136201^ 

430 

24® X 48 

3.7 

1764(x2) 

0.13620'^ 

422 

32® X 64 

4.8 

1998(x2) 

0.13632 

294 

32® X 64 

3.4 

1999(xl) 

0.13632 

289 

40® X 64 

4.2 

2028(x2) 

0.13632+ 

285 

64® X 64 

6.7 

1237(x2) 

0.13640+ 

150 

64® X 64 

3.5 

1599(x3) 


/3 = 5.40, a = 

0.060 fm, 

= 3300 MeV 


0.13640 

491 

32® X 64 

4.8 

982(x2) 

0.13647+ 

430 

32® X 64 

4.2 

1999(x2) 

0.13660 

260 

48® X 64 

3.8 

2178(x2) 


“ The number of measurements per configuration is shown in parentheses. 
1 These ensembles were generated on the QPACE systems, financed 
primarily by the SFB/TR 55, while the others were generated earlier 
within the QCDSF collaboration. 


two interpolators (33 1 to obtain an optimal interpolator Jopt = 


aJ^+P J 45 . This procedure is based on the f-dependent 2x2- 
matrix of two-point correlation functions of the interpolating 
fields J 5 and J 45 , projected onto vanishing spatial momentum. 
Solving a generalized eigenvalue problem for this matrix al¬ 
lows one to determine Jopt from the eigenvector belonging 
to the lowest energy eigenvalue. Using this interpolator in 
the correlation functions improves the signal of the ground 
state. We have also tried to apply the Additional Interpolators 
Method with a third, time-shifted interpolator ll46l l47l . but 
the results changed only marginally. Our final numbers will 
be based on the results obtained with Jopt- This differs from 
the approach of Ref. ll^ . where only the interpolator J 5 was 
utilized in the final analysis. 

To suppress statistical fluctuations, we have averaged over 
all possible values of j, k, and all possible momentum direc¬ 
tions. 


*^av 




E 

j k>j pj=±ppk=±p 


I7^: 


i l«/opf 

'4jk]4 


(57) 


where p = 27r/L. The quantities have then been fitted to 
a constant in a time interval where a plateau could be identi¬ 
fied. The choice of the fit ranges was based on the goodness of 
the correlated -values and the stability of the results upon 
reducing the fit interval. The statistical errors were evalu¬ 
ated using the Jackknife procedure combined with the binning 
method. We have observed that a binsize ribin = 4 saturates 
the statistical error, which means that the autocorrelations are 
satisfactorily taken into account. 

Our bare results are collected in Tables |V] |VI| and |VII| in 
the Appendix. In Eig. we display R~^ = (^Tyrare 
two interpolating operators J 45 and J 5 together with the cor¬ 
responding results obtained in Ref ll26l . We observe that our 
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FIG. 2. Bare results for ii^v from this work (filled symbols) and from 1261 (open symbols) for the two interpolators J45 (left panel) and J5 
(right panel). 


data are consistent with the measurements in 11 ^ . but extend 
to considerably smaller pion masses all the way down to the 
physical value. Nevertheless, in the next section we will see 
that taking into account Eq. (39 1 and using the nonperturba- 
tively computed value of C,i 2 leads to a significant shift in the 
final result. 


VI. RENORMALIZED RESULTS 


In this section we present our results for the renormal¬ 
ized quantities (43 1 ), and 

(cf. Eq. ([40)l). Eor each ensemble, the final error budget has to 
encompass the statistical errors coming from the determina¬ 
tion of the bare quantities on the lattice, the systematic uncer¬ 
tainties due to the choice of the fit range, and the errors of the 
renormalization constants. The ensuing extrapolation to the 
physical pion mass and eventually to the continuum will intro¬ 
duce further uncertainties. In order to include the errors com¬ 
ing from the renormalization constants we proceed as already 
indicated at the end of Sec.[^ Eor every fit choice in Table [II] 
we use the renormalization factors (^n, C 12 , C22 resulting from 
this fit to compute the renormalized quantities from the bare 
ratios according to Eqs. (40i and (43 1 , taking the corre¬ 
lations between and R~^ into account. The central value 


is then taken from the first fit choice, and the error due to the 
renormalization constants is determined from the differences 
with the other fit choices, as described in Sec. |IV| In the fol¬ 
lowing plots we show the central values together with their 
statistical errors, while the errors coming from the renormal¬ 
ization constants are not included, but are given in the Tables. 

We start by presenting our results for (1^)'^®. In the contin¬ 
uum limit, this quantity should be one for all pion masses. Re¬ 
sults for all ensembles are presented in Table jVIIIj In Fig. 
(l2)MS plotted for the three available lattice spacings us¬ 
ing data for tUt^L ~ 3.4 — 3.8 and ~ 260 — 294 MeV 
(or m^r 280 MeV for short; the mass dependence is rather 
weak). 



FIG. 3. ^ function of the lattice spacing a for ensembles 

with rriTrL ~ 3.4 — 3.8 and m.,r ~ 280 MeV. Only statistical errors 
are shown. 


We also show an extrapolation to the continuum limit as¬ 
suming a linear dependence on a^. We see that the result is 
consistent with unity within errors: 

= 0.9963(186)(51). (58) 

Here the first error is statistical, and the second error accounts 
for the uncertainty due to the renormalization factors, esti¬ 
mated as described at the end of Sec. jlVj It might be surpris¬ 
ing that an extrapolation linear in works so well although 
our operators are not 0(a)-improved. However, the covariant 
derivatives in the operator do not introduce 0 {a) lattice 
artifacts, at least at tree level, and the 0 {a) artifacts in O 4 
should cancel to some extent between the numerator and the 
denominator in the ratio ( |43| ). An extrapolation linear in a 
looks less stable due to the rather small range of a-values and 
yields a result which is a few percent larger. 

Note that for ^ 5 • 10“^ fm^ corresponding to /3 = 5.29, 
where most of our data are collected, we obtain, e.g., at = 
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; 

i ; 
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3 4 5 6 7 

TTlTrL 


0.3 

0.28 

0.26 



0.22 

0.2 

0.18 



0 2 4 6 8 


(10 ^ fni^) 


FIG. 4. Renormalized results af"® (upper panel) and (lower 

panel) as a function of m-nL for ensembles with /3 = 5.29 and ~ 
290 MeV. Only statistical errors are shown. 


FIG. 5. Lattice spacing dependence of (upper panel) and 
(^ 2 )MS panel) for ~ 280 MeV and m^L ~ 3.4 — 3.8. 

Only statistical errors are shown. 


294 MeV on a 32® x 64-lattice 


(n“ 0.07 fm = 0.9402(66)(54). 


(59) 


The deviation from unity is only 6%, however, it results in 
a 25 — 30% increase in the value of at the same lattice 
spacing, calculated using Eq. instead of the continuum 
relation in Eq. ([44|). 


The results for (^ 


2\MS 


and are given in Tables 


IX 


and 


[X] where the first error is statistical and the second comes 
from the uncertainty in the determination of the renormaliza¬ 
tion constants. Ideally, one would now take the infinite vol¬ 
ume limit, perform the continuum extrapolation at fixed pion 
masses and finally extrapolate to the physical mass, if it is not 
included in the range of simulated masses. Unfortunately, our 
present set of data does not allow us to perform all three ex¬ 
trapolations in a controlled way. 

We can however study the finite size effects using the data at 
P = 5.29, K = 0.13620 (m^ - 425 MeV) and k = 0.13632 
(m^ ^ 290 MeV), where we have two and three volumes, re¬ 
spectively. In Eig.j^we plot and versus for 

~ 290 MeV and see that there are indications of nonneg- 
ligible effects. In leading order chiral perturbation theory, on 
the other hand, there are no finite volume correction terms, as 
follows from the results in Ref. Il48l . 

Similarly, we use our ensembles at ^ 280 MeV and 
^ 425 MeV, where we have three and two different lattice 
spacings, respectively, to study discretization effects. Results 


for a. 


MS 


and (^ 2 )MS shown in Fig.|^ Unfortunately, wi 


with 


only three lattice spacings at hand and relatively large statis¬ 
tical errors, it is impossible to perform a reliable continuum 
extrapolation. 

According to Ref. Il48l . jjence also a^®, do not 

contain chiral logarithms, at least to one-loop order. Therefore 
we assume a linear dependence on for the extrapolation in 
the pion mass to the physical value. Since the ensemble with 
the lightest pion is already very close to the physical point, 
the chiral extrapolation is reliable. As our lattice spacings do 
not vary that much, and a proper continuum extrapolation of 
(^ 2 )MS qMS cajjjjot be attempted, we include results from 
all lattice spacings, but take into account only the data for the 
largest volume, where different volumes are available. The re¬ 
sulting extrapolations of o^® and (^ 2 ^ms jq jj^g physical pion 
mass are plotted in Fig. ^ As in these fits x^/dof is greater 
than one, we follow the procedure advocated by the Particle 
Data Group and multiply the errors by y^x^/dof. As 
before, errors coming from the renormalization constants are 
not included in the plot. We perform an extrapolation for ev¬ 
ery fit choice given in Table and compute the error of the 
final number caused by the uncertainties of the renormaliza¬ 
tion factors from the differences of the extrapolated results as 
indicated at the end of Sec. m 

From this procedure we find our final results 

(^2^MS ^q 2361(4i)(39)^ 

of® = 0.1364(154)(145) (60) 

at the scale p, = 2 GeV. They can be compared with the 
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FIG. 6. Extrapolation to the physical pion mass for al*® (left panel) and panel). The open triangle represents the extrapolated 

value. Only statistical errors are shown. 


earlier lattice calculations 

(^2^MS ^ 0.269(39), ap = 0.201(114), ED 

(^2)MS ^ 0.28(1)(2), af®= 0.233(29)(58), E?) 

( 61 ) 

where, for ETII . we have quoted the result for (^ 2 ^ms qjj jjjgjj. 
larger lattice and used the continuum relation in Eq. ( |44) l to 
calculate the corresponding value of the second Gegenbauer 
moment 

It should, however, be kept in mind that all these numbers 
were obtained on lattices with lattice spacings between 0.06 
and 0.08 fm. The investigation of discretization effects for 
(^ 2 )MS ^MS remain a challenge for future studies. 


VII. CONCLUSIONS AND OUTLOOK 

We have presented the most accurate, up to now, lattice 
determination of the second moment of the pion distribution 
amplitude using two flavors of dynamical (clover) fermions 
on lattices of different volumes and pion masses down to al¬ 
most the physical value. So the chiral extrapolation per se 
does not seem to be an issue. Also the omission of strange 
quarks should not be of great importance. However, the sta¬ 
tistical fluctuations of the lattice matrix elements of operators 
with derivatives are large for small pion masses and require 
averaging over a large number of configurations in order to 
obtain phenomenologically relevant precision. We found that 
the signal can be somewhat improved by using the variational 
method with the two interpolators corresponding to the pseu¬ 
doscalar and axial-vector currents. 

The main difference of this work from the previous stud¬ 
ies EhllZTl is the nonperturbative evaluation of the full 2x2 
mixing matrix of the operators with two derivatives. In the 
framework of Ref. E6l the nonperturbative mixing coefficient 
turns out to be of the same sign but up to one order of magni¬ 
tude larger than the same coefficient computed perturbatively. 


This observation underlines the necessity of nonperturbative 
renormalization, at least at the presently reachable /3-values. 

Still, some uncertainty in the renormalization factors re¬ 
mains. It is dominated by the uncertainty in the conversion 
factors connecting the RT-SMOM scheme to the MS scheme, 
which are calculated in continuum perturbation theory and are 
known to two-loop accuracy ESI Eh). A three-loop calcula¬ 
tion is, therefore, needed in order to further reduce the renor¬ 
malization uncertainty and would be extremely welcome. 

In our work we have also emphasized the importance of 
using the corrected relation Eq. (451 between ^ms 

for finite lattice spacing, instead of the continuum relation in 
Eq. ([g, due to discretization errors in derivatives that lead to 
a violation of the product rule. This effect is studied in detail. 

From our data we cannot exclude significant discretiza¬ 
tion effects in (^ 2 ^ms a^®, but a quantitative study re¬ 

quires simulations at smaller lattice spacings of the order of 
a ~ 0.04 fm, which are presently not available to us. Such lat¬ 
tices wifi be generated in the future within the CLS effort ifSOl . 
This will be a major step towards the calculation of the second 
moment of the pion DA with fully controllable accuracy. As a 
final remark, we note that the somewhat smaller value of 
obtained in this work seems to be favored by the phenomeno¬ 
logical studies of form factors in the framework of light-cone 
sum rules, see, e.g.. Refs. QiiiiiHia. 
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collaboration for giving us permission to use some of their Appendix A: Bare and renormalized results by the ensemble 
gauge field configurations. 

The following Tables summarize the results obtained for 
each gauge field ensemble separately. 


TABLE V. Bare results for using J 5 as interpolator. 


p 

Hi 

Size 

Fit range 

p — 

-^av 

xVdof 

Fit range 

p+ 

-^av 

xVdof 

5.20 

0.13596 

32® X 64 

3-14 

0.1674(36) 

0.67 

9-19 

0.6013(46) 

6.67 

5.29 

0.13620 

24® X 48 

8-12 

0.161(25) 

0.47 

7-12 

0.5792(97) 

4.44 

5.29 

0.13620 

32® X 64 

3-14 

0.1668(30) 

0.73 

11 - 20 

0.6187(52) 

7.73 

5.29 

0.13632 

32® X 64 

3-17 

0.1705(78) 

0.72 

9-15 

0.602(11) 

7.76 

5.29 

0.13632 

40® X 64 

3-18 

0.1756(33) 

1.51 

10-25 

0.6213(36) 

5.38 

5.29 

0.13632 

64® X 64 

7-15 

0.1694(37) 

0.82 

15-26 

0.6343(22) 

6.67 

5.29 

0.13640 

64® X 64 

5-20 

0.1627(56) 

0.76 

17-25 

0.6421(61) 

6.00 

5.40 

0.13640 

32® X 64 

3-15 

0.1679(38) 

0.59 

14 - 25 

0.654(14) 

3.83 

5.40 

0.13647 

32® X 64 

3-15 

0.1653(35) 

1.04 

15 - 22 

0.657(21) 

2.78 

5.40 

0.13660 

48® X 64 

3-15 

0.1681(32) 

0.81 

15-25 

0.6467(57) 

4.22 


TABLE VI. Bare results for using J 45 as interpolator. 

P 

Hi 

Size 

Fit range 

r>— 

-^av 

xVdof 

Fit range 

p+ 

-^av 

xVdof 

5.20 

0.13596 

32® 

X 

64 

10-16 

0.1859(91) 

1.57 

13-19 

0.6354(72) 

1.29 

5.29 

0.13620 

24® 

X 

48 

7-13 

0.1845(83) 

1.08 

10-15 

0.680(12) 

0.92 

5.29 

0.13620 

32® 

X 

64 

10-15 

0.1963(60) 

0.13 

18 - 24 

0.617(10) 

1.86 

5.29 

0.13632 

32® 

X 

64 

9-15 

0.155(14) 

0.25 

12 - 20 

0.660(13) 

0.67 

5.29 

0.13632 

40® 

X 

64 

8-15 

0.1976(49) 

0.83 

17-24 

0.6441(74) 

1.19 

5.29 

0.13632 

64® 

X 

64 

10-25 

0.1839(39) 

1.59 

16-30 

0.6394(23) 

1.64 

5.29 

0.13640 

64® 

X 

64 

10-19 

0.2015(97) 

0.77 

20-30 

0.6321(67) 

0.37 

5.40 

0.13640 

32® 

X 

64 

7-15 

0.1931(42) 

0.41 

16-25 

0.682(13) 

0.77 

5.40 

0.13647 

32® 

X 

64 

3-13 

0.1980(17) 

0.59 

17-22 

0.682(17) 

0.22 

5.40 

0.13660 

48® 

X 

64 

14 - 20 

0.1823(93) 

0.60 

19-29 

0.6640(73) 

0.52 



TABLE VII. Bare results for 

. using the variational method with the intrpolators J 45 , J 5 . 


P 

K 

Size 

Eit range 

-^av 

xVdof 

Fit range 

p+ 

-^av 

xVdof 

5.20 

0.13596 

32® 

X 

64 

3-16 

0.1813(27) 

0.63 

10-19 

0.6142(46) 

0.52 

5.29 

0.13620 

24® 

X 

48 

3-13 

0.1660(52) 

1.01 

5-13 

0.6039(54) 

0.38 

5.29 

0.13620 

32® 

X 

64 

4-16 

0.1775(32) 

0.52 

9-16 

0.6303(35) 

0.41 

5.29 

0.13632 

32® 

X 

64 

6-16 

0.1710(120) 

0.63 

5-16 

0.6289(44) 

0.35 

5.29 

0.13632 

40® 

X 

64 

2-23 

0.1838(24) 

1.52 

14 - 24 

0.6226(56) 

0.40 

5.29 

0.13632 

64® 

X 

64 

2-22 

0.1761(21) 

0.85 

8-25 

0.6353(14) 

0.93 

5.29 

0.13640 

64® 

X 

64 

2-20 

0.1790(39) 

0.78 

10-20 

0.6350(30) 

1.35 

5.40 

0.13640 

32® 

X 

64 

2-14 

0.1773(27) 

0.55 

13-20 

0.657(11) 

0.45 

5.40 

0.13647 

32® 

X 

64 

2-16 

0.1742(22) 

1.03 

16 - 22 

0.662(25) 

0.26 

5.40 

0.13660 

48® 

X 

64 

2-16 

0.1794(24) 

0.80 

15-25 

0.6534(53) 

0.30 
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TABLE VIII. Results for = 2 GeV) using the variational 

method with the interpolators J 45 , J 5 . The first error corresponds 
to the statistical fluctuations, and the second to the contribution from 
the uncertainty in the determination of the renormalization constants. 

P 

K 

Size 

{l®)®®(/r = 2GeV) 

5.20 

0.13596 

32® X 64 

0.9298(70) (56) 

5.29 

0.13620 

24® X 48 

0.9028(81)(52) 

5.29 

0.13620 

32® X 64 

0.9422(53) (55) 

5.29 

0.13632 

32® X 64 

0.9402(66) (54) 

5.29 

0.13632 

40® X 64 

0.9308(84) (54) 

5.29 

0.13632 

64® X 64 

0.9498(20) (55) 

5.29 

0.13640 

64® X 64 

0.9494(44) (55) 

5.40 

0.13640 

32® X 64 

0.9690(159)(51) 

5.40 

0.13647 

32® X 64 

0.9757(371)(51) 

5.40 

0.13660 

48® X 64 

0.9632(79) (50) 


TABLE IX. Results for (C^)^®(/i = 2 GeV) using the variational 
method with the interpolators J 45 , J 5 . The first error corresponds 
to the statistical fluctuations, and the second to the contribution from 
the uncertainty in the determination of the renormalization constants. 

P 

K, 

Size 

(e®)"'®(/r = 2 GeV) 

5.20 

0.13596 

32® X 64 

0.2427(53)(28) 

5.29 

0.13620 

24® X 48 

0.2147(103)(42) 

5.29 

0.13620 

32® X 64 

0.2325(63)(42) 

5.29 

0.13632 

32® X 64 

0.2199(240)(45) 

5.29 

0.13632 

40® X 64 

0.2467(49)(38) 

5.29 

0.13632 

64® X 64 

0.2289(42)(44) 

5.29 

0.13640 

64® X 64 

0.2348(78)(42) 

5.40 

0.13640 

32® X 64 

0.2284(58)(49) 

5.40 

0.13647 

32® X 64 

0.2212(66)(51) 

5.40 

0.13660 

48® X 64 

0.2334(50)(48) 


TABLE X. Results for = 2 GeV) using the variational 

method with the interpolators J 45 , J 5 . The first error corresponds 
to the statistical fluctuations, and the second to the contribution from 
the uncertainty in the determination of the renormalization constants. 

P 

K 

Size 

= 2 GeV) 

5.20 

0.13596 

32® X 64 

0.1654(161)(113) 

5.29 

0.13620 

24® X 48 

0.0996(304) (152) 

5.29 

0.13620 

32® X 64 

0.1286(188)(153) 

5.29 

0.13632 

32® X 64 

0.0930(700)(161) 

5.29 

0.13632 

40® X 64 

0.1767(160)(141) 

5.29 

0.13632 

64® X 64 

0.1136(124)(158) 

5.29 

0.13640 

64® X 64 

0.1310(230)(154) 

5.40 

0.13640 

32® X 64 

0.1010(214)(169) 

5.40 

0.13647 

32® X 64 

0.0760(380) (176) 

5.40 

0.13660 

48® X 64 

0.1188(159)(164) 
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